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Abstract 

«j . We consider N identical oscillators coupled to a single environment and show that the conditions for 

' the existence of decoherence free subspaces are degeneracy of the oscillator frequencies and separability 

of the coupling with the environment. A formal exact equation for the evolution in the case of two 
£SJ ■ oscillators is found and the decoherence free subspace is explicitly determined. A full analytical solu- 

tion for any initial condition and general parameters (frequencies and dissipation constants) is given in 
the markovian approximation and zero temperature. We find that slight relaxation of degeneracy and 
separability conditions leads to the appearence of two components in the dynamical evolution with very 
different decoherence times. The ratio between the characteristic time of the weak and strong decoherent 
components is given by twd/tsd ~ k/5k, where 5k is a measure of the nonseparability of the coupling 
■ to the environment and k is the mean decay constant of the oscillators. 

3 ■ 



1 Introduction 



The very same mechanism responsible for the potential improvements on computation speed using quantum 
mechanics, is the one which greatly hinders immediante technical implementation. Entanglement between 
different subsystems is essential for the production of the states used in information processing; at the 
same time it prevents these qubits to be completely isolated from its environment, producing undesired 
entanglement with the environmental degrees of freedom. The deleterious effect of this coupling is usually 
called decoherence. Therefore much effort has been devoted to finding ways around decoherence in quantum 
computation, such as error correcting codes (^], dynamical decoupling || and computation in decoherence 
free subspaces [f|. Experimental observations of decoherence free evolution have been reported. (5[ ^| Many 
physical implementations have been proposed including cavity QED[Q, ions traps nuclear magnetic 
resonance Q and semiconductor quantum dotsjio||. From the theoretical point of view, recent work has 
been mainly on proving the existence of DF subspaces. Much less work, however has been devoted to the 
dynamics. In general those models are highly idealized, and it is of interest at least to relax some of the 
stringent conditions, e.g. degeneracy of the free modes, and calculate such effects on the time evolution of 
the system. This could be considered as an important step towards realistic implementations of the models. 

In the present work we consider the case of N independent oscillators linearly coupled to a single envi- 
ronment, and show that the existence of strict decoherence free subspaces can only be obtained under the 
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following two conditions: degeneracy of the oscillators, and separability of the coupling with the environ- 
ment. Both the exact form of the spectral density and the temperature of the environment are immaterial 
in what concerns the existence of DFS. A formal equation for the dynamics of the reduced density matrix 
of the oscillators is derived which resembles the corresponding master equation, but with time dependent 
coefficients. For two independent oscillators we solve the dynamics of the reduced density and construct the 
decoherence free subspace. Moreover under Markov approximation, and zero temperature, we are able to 
give an analytical expression for the time evolution of the system of interest, in the general case (i.e. non 
degenerate modes and nonseparable coupling) . We study the effect of a slight breaking of the degeneracy and 
separability and show that the dynamics will in general produce a long lived and a short live components. 
The time scale for the duration of these components is derived in terms of the appropriate parameters. 



2 The model 

Let us consider a collection of N identical harmonic oscillators linearly coupled to a single environment, 
modelled as a large set of harmonic oscillators, 
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If the coupling of the N oscillators to the environment does not depend on the particular oscillator, apart 
from an overall factor, the coupling constants gik can be written as a product GiD^. In this case it is 
immediate to verify that there is a single collective mode which will remain coupled to the environment, 
while the other N-l collective modes form a decoherent free (infinite) subspace. j| The collective mode can 
be easily visualized if we rewrite the hamiltonian as follows 

H = W ^a 1 ta i + ^^.bJ c b k + ^ (dI(Y,G*^)W + D^G^bn , 

i=l fe k \ i i I 

where it is clear that the collective mode coupled to the environment has the following creation and annihi- 
lation operators 

it _ J2i G*a\ ~ _ J^i GjM 

This can be viewed as a N-dimensional rotation in the space of the original creation and annihilation oper- 
ators, which maps the set ai onto the set Ai in terms of which the hamiltonian is rewritten as 

N 

H = u ^2 k\ Ai + Y ^feb k b k + Ai Y cfeb k + A x ^ cj6£, 

i fe fe fe 

with Cfc = Yli \Gi\ 2 D^. The rotation leading to the equation above defines the relation between the new and 
old operators. This rotation is not unique except in the case of two identical operators. We show next that in 
contrast with other commonly studied models, in the present one we are able to introduce both decoherence 
and dissipation, for arbitrary environmental spectral densities (ohmic, subohmic, etc) and temperatures. A 
long but straighforward procedure leads to the following exact generalized master equation 



dp _ 1 
dt ih 



(X + e)(2A 1 -Al-A\A 1 ---A\A 1 )p 



H ,P 

e(2A[ ■ Ai - AiA[- - 'AjA^ 
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where 



A' 



H = fiw^AtAi + ^w + ^AiAi. 

i=2 

The real functions A, S, e are implicitly defined in terms of the auxiliary function r)(t) 

r](t) = cxp (- X(t')dt' -iwt-ij 5(t')dt'^j 



which satisfies the integrodifferential equation 



77 + iujr] + 



[ dTy2\ck\ 2 e iUk(t - T) ri(T)=0, 
Jo u 



(2) 



(3) 



subject to the initial condition 77(0) = 1. Moreover we have 



Ck\ 2 n k (/3) 

\v(t)\ 2 



d Te -*»k(t-T) v ^ 



where rik(/3) is the mean excitation number for the k-th mode of the environment at inverse temperature 
[3 = 1/ksT. If the usual Born-Markov approximations hold, then S(t) = 0, A(t) = ^fci, and e = 
where fcj characterize the markovian evolution when the only the i-th original oscillator is coupled to the 
bath, and n is the environment mean number of thermal excitations. In this case the master equation 
simplifies to 



dp 
~dl 



ico 



[AjAi, p] + k *(™ + !)( 2A i • A i - A i A i- - -M^i)/ 3 

t=l i=l 
N 

+ hn{2k\ ■ Ai - AiAi- - -AxA^p. 



(4) 



Note that the collective mode which couples to the environmental degrees of freedom decoheres much faster 
than any of the individual oscillators. 

2.1 Two Identical Oscillators 

For the sake of definiteness we restrict ourselves to the case of two identical oscillators. Let us see that in 
this case the collective mode Ai (A2) which couples (decouples) to the environment is given by 



A(t) 



cos(6»)4 t} + sin(0)ay 



(t) 



(t) 



sin(6')4 t) + cos(0)a£ 



s(t) 



with tan(0) = G2/G1, which reduces to tan(0) = \Jk%jk\ when the markovian approximation is valid. 
Let us see that the transformation above can be implemented as using an appsropriate rotation operator. 
The generalized master equation (|l|) for two identical oscillator can be exactly solved using the evolution 
superoperator U (t) 



where the coefficients v(i), xit) and z(t) can be given in terms of the functions rj(t), cq. (0), and Af(t), 



(5) 



Af(t) = / dre(r) 
Jo 



77(1-) 
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as follows 



w 1+Af(t)' w y/l+Af(t) 1+Af(t) 



In the markovian limit the preceding formulas reduce to 



1 e (-iui-k l ~k 2 )t (n + 1)(1 — e~ 2 ( kl+k2 ' )t ) 

x = In — , 2 — 



1 + n{l - e - 2 ( fc i+ fc =)*) ' v / 1 + fl (i_ e -2(fc 1 +fe 2 )t)' l + n(l-e- 2 ( fe i+ fe 2)t)) ' 

In eq. (|J) the first and last terms of the rhs correspond to the rotation which leads to the coupled and 
uncoupled collective modes. 

Since the second collective mode is effectively decoupled from the environment, any density operator 
in the Hilbert space of this mode, times the asymptotic density operator of the coupled collective mode, 
provided it exists, will experience a unitary evolution. For simplicity we will further on restrict ourselves to 
the zero temperature case. Any density operator of the form 

3 = V ^n^(-al sin(0) + a* cos(6»)) n |0, 0)(0, 0|(-ai sin(0) + a 2 cos(6»)) m 
^-^ n\m\ 



^ . V^!^(-sm60™ +m_ni ~ mi (cos6>) ni+mi 
nn, mi ^ y/(n-n 1 )\n 1 \(m-m 1 )\m 1 \ 
x |m, n — ni)(mi, m — mi\, 

will be protected against dissipation and decoherence. In fact, applying the evolution superoperator to an 
initial density matrix of this form, we obtain 



Pit) 



E 



-iujt(n-m) * 



V n!m!(— sin#) 



n-\-m—n\ —m 



1 (cos I 



\n 1 +m 1 



n,m,ni ,mi 



y{n — rti)!rti!(m — mjlmi! 
x \m, n — ni)(mi, m — m\\, 



Now, we use the evolution superoperator on the initial operator density 

p(0) = (cos(a)|l,0) +sin(a)e^|0, l))(cos(a)(l,0| +sin(a)e-^(0, I|), 



(6) 



to obtain the time dependent evolution, which can be given in closed form and corresponds to the limit 
where the frequencies are degenerate and coupling is separable, of the formulas given in the next section. 
For now we show the asymptotic limit, 

p(t * oo) = P\1>)W + (1 - P)|0, 0) (0, 0| 

where the probability to go to the state \ip), 



is given by 



P 



/fc^|i,o) -Vh\o,i) 



fk~2 cos(a) — vfci e sin(a) 



y/ki + k 2 

and the probability to go to the joint ground state is 1 — P. 
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Observe that varying a and <j> we can go from total preservation to total leakage. For example, if we set 
tan(a) = ^/k\/k2, and <f> = then the full state will leak to the ground state |0,0). On the other hand, if 
we set tan(a) = —■s/ki/ki, and <f> — then the initial state will be exactly equal to |?/>), and will persist at 
all times with probability 1. All other combinations will go to the asymptotic state \ip) with probability P 
and to the ground state 1 0, 0) with probability 1 — P. The asymptotic fidelity, F(oo), which is the overlap 
between the initial and final density matrices, is given by 

(v / Ai2Cos(a) — \fk[e 1 ^ sin(a))(v / fc2 cos(a) — \fk[e~ 1 ^ sin(a)) 
~ ki +k 2 



3 Effects of more Realistic Modelling 

We remark that the results above were obtained under a number of assumptions, which will be relaxed below. 
Notice that the use of the rotating wave approximation (RWA) is not essential in obtaining the decoupled 
mode: any interaction linear in the field operators would be as good (provided the other assumptions hold) . 
Had we chosen an interaction linear in the identical oscillators but nonlinear on the environmental operators 
we would have obtained also a decoupled collective mode. In these cases, however, the complication would 
be only of technical nature leading to (much) more complex dynamics. 

Another important hypothesis to obtain DFS is that of identical frequencies of the original main oscilla- 
tors. Of course, any interaction between them would destroy the symmetry upon which the existence of DFS 
rests. On the other hand, we have assumed that the oscillator-environment coupling satisfies = GiD^, 
which amounts to a separable coupling. It is not an easy task to find realizations of such interactions in 
nature given its nonlocal character. However, it might be a good approximation in special circumstances, as 
e.g. optical cavities. A particular consequence of the separability hypothesis can be seen writing the master 
equation, in the zero temperature limit, in terms of the original oscillators (with different frequencies for 
generality) 

Cq = (— iu>± — k\) + -a|a (iu>\ — k\) + 2fcxav*4 

+a 2 a 2 - (-kj 2 - k 2 ) + -a 2 a 2 (iuj 2 - k 2 ) + 2fc 2 a 2 -a 2 

2ai-a 2 — a|a 2 - — -a|a 2 ) + k 3 (2a 2 -a| — aa 2 - — ai-a 2 ) . (7) 



+k 3 (5 



Notice that we have introduced a new quantity k 3 , which is the environment induced coupling. When the 
separability condition holds the three dissipation constants, fci,fc 2 ,fc3 satisfy the relationship fcf = fci& 2 . 
Since these constants are given by the expressions 

i i2 I \2 k% » . . k 3 

\9ik\ ~ — , \92k\ ~ — , 3ifcff 2 fcM~— (8) 

TT TT TT 

Schwarz inequality implies immediately that \k 3 \ 2 < k\k 2 in the general (nonseparable) case. In what follows 
we present the general solution to eq. (Q) and study, in particular, the dynamics for parameters which do 
not quite satisfy the condition for DFS, but both the frequency degeneracy and coupling separability are 
almost fulfilled. 

The evolution superoperator for eq.(Q) is given by 1 12 



p(t) — e J1 ^ ai ' a i e : ' 2( - t - )a2 ' a 2 g z (*) a 2- a i e z " '(t)ai-a2 e 9(*) a l a 2' e 9*(*)- a i a 2 e m 2( t ) a 2 a 2- e m 2(*)- a 2 a 2 Q 

0g"»i (*) a l a i- e m i (t)- a l a l e ?(*) a i a 2- R q* (*)ai -] ea 2p (9) 
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where 



R = 



~2 + 2 • 

q(t) = 2k 3 (l-e rt ) (A + e r 



c = k 2 — k\ + i (uJ2 — u\) , r = <J c? + 4fc|, A± = c ± r 



for r ^ 



o m l(*) 



-i?, t 



2r 



-e~V (A+e rt 



j 2 (i) = (l + | g (i)| 5 



,m 2 (t) 



- 1 



ix(t)= e- mi W + g(t) 2 e-"* a W + (| g (t)| 2 ) 

*(t) = - (? (t) e -(™K*)+™2(*)) - g*( t ) (1 + |g(t)| 2 ) 



e ™ 2 (t) _ e -2Bt e -m 1 (i) 



,7712 (t) 



,77(2 (t) 



Note that the neither the separability nor the degeneracy conditions have been used so far in obtaining the 
solution. For the sake of comparison we use the same initial condition of the precedent section (eq. @). In 
the general case its time evolution is given by 

Pit) = p(t)\m)(m\ + (i - p(i))io,o><o,oi 

where the probability to go to the state \ip(t)), 

(cos(0)M_(i) +sm(0)e i<#, QO))|l,O) + (sin(6»)e^Af + (t) + cos(6»)Q(i))|0, 1) 



is given by 



W)) = 



P(t) = ||cos(0)M_(f) +sm(6)e i ' l> Q(t)\ 



3in(0)e^M + (t) + cos(0)Q(i)|| , 



and the probability to go to the joint ground state is 1 — P(t). The functions M±(t) and Q(t) are given by 



M ± (t) 



-m 



-rt/2 nzc C t/2 < 



2 (lT-) + e rt ^(l±-) , Q{t) 



(e~ rt 



/2 _ e rt/2 



Now, we assume slight deviations from degeneracy and separability, that is, u>% — to — 8u>, ll>2 — uj + 5u, 
^3 = \fk\k2 — 5k, with 5u> <C u), 5k -C s/kxkz. Then, the state \ip{t)} can be approximated as 

\m) = 

where Q , £j do not have any temporal dependence and are given by 



ut x / Ci|i,Q) +6IQ.1) r - ffcl+ ^ t GMH»i) e -^ 



(fc^i-j =p i&j) cos(a) ± vfcife sin(a) 



3*^ 



(A^ 2 -j ± iSuj)e 1 ^ sin(a) ± \Jk\ki cos (a) 



(10) 



(11) 



In the general case it is not possible to find initial conditions which are completely decoherence free. Nor it 
is possible to find to orthogonal subspaces with very different characters in what decoherence is concerned. 
However, we can choose the initial condition as to have a minimal component either in a strong decoherence 
(SD) or in a weak decoherence (WD) subspaces, by choosing, e.g. 



tania 



)/SD\ 



(l) 



5lu 



6 



where k is some average dissipation constant. The corresponding subspaces, apart from a phase, can be 
written as 



\Jk\ + k 2 

If Su> <C k\, k 2 the phase can be ignored. Moreover, if k\ — k 2 then we have k = k\ = k 2 and the phase can 
be unambigously determined. The weak decoherence wavefunction defines a space which is robust against 
decoherence with typical decoherence time scale of the order 

fci + k 2 

twd = 



which is to be compared with the strong decoherence time associated to \^Psd) given by 

1 

T.SD 



ki + k 2 



If ki and k 2 are of approximately the same order (fc), then the weak decoherence time is about k/Sk times 
greater than the strong decoherence time. In all other cases (different from Sui = or ki = k 2 one always 
has, for any initial condition (of the form @), a term which decays slowly and one which decays rapidly 
comparatively. However, the forms given above may help one to choose the initial condition in such a way 
that the time evolution is robust against decoherence (it should be close to form ( |l2| ) with k — (k\ + k 2 ) /2). 

Had we chosen other Fock subspaces (for example those of m total excitations), we would have obtained 
also an analytical solution showing both slowly and fast decaying components. Of course, the expressions 
get more involved. 



4 Concluding Remarks 

We have studied decoherence and dissipation free subspaces for systems involving N identical harmonic 
oscillators linearly coupled to a single environment. The conditions for the existence of such decoherent free 
subspaces are too stringent, and hardly implementable in practice. Therefore, we studied the dynamics of 
such systems when these conditions are slightly broken. This may serve as a guide for preparation of initial 
conditions wich are fairly robust against decoherence. To the extent that this model is adequate for the 
description of modes in a QED cavity, it would be possible to use it to generate good quantum memories in 
the presence of dissipation and decoherence. 
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